*MAIN ANALYSIS DOFILE

clear
eststo clear

*cd "~/Dropbox/Donor Adjustment/Replication/"
use "CombinedData.dta", clear

order WDI_Name id year
xtset id year

*** DATA MANIPULATION

gen log_SO2 = log(SO2+1)
label var log_SO2 "log (SO2)"
replace CO2=0 if CO2 < 0
gen log_CO2 = log(CO2+1)
label var log_CO2 "log (CO2)"
gen log_BOD = log(BOD)
label var log_BOD "log (BOD)"
gen log_gdp = log(gdp)
gen log_PM10 = log(pm10)
label var log_PM10 "log (PM10)"
label var log_gdp "log (GDP)" 
gen log_gdp_pc = log(gdp_pc)
label var log_gdp_pc "log (GDP p.c.)"
gen pop = gdp/gdp_pc
gen log_pop = log(pop)
label var log_pop "log (Population)"
label var polity2 "Polity"
gen time = year-1960
gen time2 = time^2
gen time3 = time^3
gen wto_polity2 = wto * polity2
encode(region), gen(region_id)

gen log_deforest=log(deforest+1)

gen envmin=env_ministry_real
by id (year), sort: replace envmin =. if envmin[_n-1]== 1 
by id (year), sort: replace envmin =. if envmin== 1 & envmin[_n-1]==.

gen FIT1=FIT
by id (year), sort: replace FIT1 =. if FIT1[_n-1]== 1 
by id (year), sort: replace FIT1 =. if FIT1== 1 & FIT1[_n-1]==.
replace FIT1=1 if id==72 & year==1990

*** DYNAMIC MODEL FOR POLLUTION

gen L1log_SO2=L1.log_SO2
label var L1log_SO2 "log(SO2) (t-1)" 
gen L1log_CO2=L1.log_CO2
label var L1log_CO2 "log(CO2) (t-1)"
gen L1log_BOD=L1.log_BOD
label var L1log_BOD "log(BOD) (t-1)"
gen L1log_PM10=L1.log_PM10
label var L1log_PM10 "log(PM10) (t-1)"

gen L1gatt_wto=L1.gatt_wto
label var L1gatt_wto "GATT/WTO Membership (t-1)"
gen L1log_gdp_pc=L1.log_gdp_pc
label var L1log_gdp_pc "log(GDP p.c) (t-1)"
gen L1log_pop=L1.log_pop
label var L1log_pop "log(population) (t-1)"
gen L1life_expectancy=L1.life_expectancy
label var L1life_expectancy "Life expectancy (t-1)"
label var life_expectancy "Life expectancy"
gen L1polity2=L1.polity2
label var L1polity2 "Polity (t-1)"
gen L1log_environment=L1.log_environment
label var L1log_environment "log(Environmental Aid) (t-1)"
gen L1log_neutral=L1.log_neutral
label var L1log_neutral "log(Neutral Aid) (t-1)"
gen L1log_dirty=L1.log_dirty
label var L1log_dirty "log(Dirty Aid) (t-1)"
label var urban_pop_per "Urban Population $\%$" 
gen L1urban_pop_per=L1.urban_pop_per
label var L1urban_pop_per "Urban Population $\%$ (t-1)"
gen L1electricity_coal_production=L1.electricity_coal_production 
label var L1electricity_coal_production "Electricity from Coal $\%$ (t-1)" 
label var electricity_coal_production "Electricity from Coal"
gen L1electricity_production_renew=L1.electricity_production_renew
label var L1electricity_production_renew "Electricity from Renewable Energy $\%$ (t-1)" 
label var electricity_production_renew "Electricity from Renewables"
gen L1wto=L1.wto
label var L1wto "WTO Membership (t-1)"
label var wto "WTO Membership"

gen L1gatt_wto_informal=L1.gatt_wto_informal
label var L1gatt_wto_informal "GATT/WTO Informal Standing (t-1)"
label var gatt_wto_informal "GATT/WTO Informal Standing"

egen PTA = rowtotal(USPTA JapanPTA EFTAPTA EUPTA)

gen post_tokyo = 0 
replace post_tokyo = 1 if year > 1979
label var post_tokyo "Post-Tokyo Round"

gen post_uruguay = 0 
replace post_uruguay = 1 if year > 1994
label var post_uruguay "Post-Uruguay Round"

gen L1post_tokyo_gatt_wto = L1gatt_wto * L1.post_tokyo
label var L1post_tokyo_gatt_wto "GATT/WTO Membership * Post-Tokyo (t-1)"

gen post_tokyo_gatt_wto = gatt_wto * post_tokyo
label var post_tokyo_gatt_wto "GATT/WTO Membership * Post-Tokyo"

gen L1post_uruguay_gatt_wto = L1gatt_wto * L1.post_uruguay
label var L1post_uruguay_gatt_wto "GATT/WTO Memberhsip * Post-Uruguay (t-1)"

gen post_uruguay_gatt_wto = gatt_wto * post_uruguay
label var post_uruguay_gatt_wto "GATT/WTO Memberhsip * Post-Uruguay"

gen L1auto_emissions_reg=L1.auto_emissions_reg
label var L1auto_emissions_reg "Auto Emission Standards (t-1)"


*cd "~/Dropbox/Donor Adjustment/Replication/check"

********************************************
*****Domestic Political Institutions********
********************************************

*Table A37: Difference in Difference: Estimating the Impact of GATT/WTO on SO2 by Regime Type
set more 1
eststo clear
eststo: xi: xtreg log_SO2 gatt_wto i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto log_gdp_pc log_pop  i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_SO2 gatt_wto i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_SO2 gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_SO2 wto i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_SO2 wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_SO2 wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust

eststo: xi: xtreg log_SO2 wto i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_SO2 wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_SO2 wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
*esttab using "SO2Regime(DID).tex", booktabs label title(Difference in Difference: Estimating the Impact of GATT/WTO on $SO2$ by Regime Type \label{so2(DID)Demo}) keep(gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A38: Dynamic Model: Estimating the Impact of GATT/WTO on SO2 by Regime Type
eststo clear
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto i.year i.region  if polity2 > 6 & polity2 != ., robust 
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != ., robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != ., robust

eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto i.year i.region  if polity2 <= 6 & polity2 !=., robust 
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != ., robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != ., robust

eststo: xi: xtreg log_SO2 L1log_SO2 L1wto i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust

eststo: xi: xtreg log_SO2 L1log_SO2 L1wto i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_SO2 L1log_SO2 L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
*esttab using "SO2Regime(Dynamic).tex", booktabs label title(Dynamic Model: Estimating the Impact of GATT/WTO on $SO2$ by Regime Type\label{so2(Dynamic)Demo}) keep (L1log_SO2  L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_SO2 )  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A39: Difference in Difference: Estimating the Impact of GATT/WTO on CO2 by Regime Type
eststo clear
eststo: xi: xtreg log_CO2 gatt_wto i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_CO2 gatt_wto i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_CO2 gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_CO2 wto i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_CO2 wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_CO2 wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust

eststo: xi: xtreg log_CO2 wto i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_CO2 wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_CO2 wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
*esttab using "CO2Regime(DID).tex", booktabs label title(Difference in Difference: Estimating the Impact of GATT/WTO on $CO2$ by Regime Type \label{CO2(DID)Demo}) keep(gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A40: Dynamic Model: Estimating the Impact of GATT/WTO on CO2 by Regime Type
eststo clear
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto i.year i.region  if polity2 > 6 & polity2 != ., robust 
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != ., robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != ., robust

eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto i.year i.region  if polity2 <= 6 & polity2 !=., robust 
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != ., robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != ., robust

eststo: xi: xtreg log_CO2 L1log_CO2 L1wto i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust

eststo: xi: xtreg log_CO2 L1log_CO2 L1wto i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_CO2 L1log_CO2 L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
*esttab using "CO2Regime(Dynamic).tex", booktabs label title(Dynamic Model: Estimating the Impact of GATT/WTO on $CO2$ by Regime Type\label{CO2(Dynamic)Demo}) keep (L1log_CO2  L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_CO2 )  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A41: Difference in Difference: Estimating the Impact of GATT/WTO on BOD by Regime Type
eststo clear
eststo: xi: xtreg log_BOD gatt_wto i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_BOD gatt_wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_BOD gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_BOD gatt_wto i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_BOD gatt_wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_BOD gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_BOD wto i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_BOD wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_BOD wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust

eststo: xi: xtreg log_BOD wto i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_BOD wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_BOD wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
*esttab using "BODRegime(DID).tex", booktabs label title(Difference in Difference: Estimating the Impact of GATT/WTO on $BOD$ by Regime Type \label{BOD(DID)Demo}) keep(gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A42: Dynamic Model: Estimating the Impact of GATT/WTO on BOD by Regime Type
eststo clear
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto i.year i.region  if polity2 > 6 & polity2 != ., robust 
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != ., robust
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != ., robust

eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto i.year i.region  if polity2 <= 6 & polity2 !=., robust 
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != ., robust
eststo: xi: xtreg log_BOD L1log_BOD L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != ., robust

eststo: xi: xtreg log_BOD L1log_BOD L1wto i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust

eststo: xi: xtreg log_BOD L1log_BOD L1wto i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_BOD L1log_BOD L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
*esttab using "BODRegime(Dynamic).tex", booktabs label title(Dynamic Model: Estimating the Impact of GATT/WTO on $BOD$ by Regime Type\label{BOD(Dynamic)Demo}) keep (L1log_BOD  L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_BOD )  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A43: Difference in Difference: Estimating the Impact of GATT/WTO on PM10 by Regime Type
eststo clear
eststo: xi: xtreg log_PM10 gatt_wto i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_PM10 gatt_wto i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg log_PM10 gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != .,  fe robust

eststo: xi: xtreg log_PM10 wto i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_PM10 wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_PM10 wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust

eststo: xi: xtreg log_PM10 wto i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_PM10 wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg log_PM10 wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
*esttab using "PM10Regime(DID).tex", booktabs label title(Difference in Difference: Estimating the Impact of GATT/WTO on $PM10$ by Regime Type \label{PM10(DID)Demo}) keep(gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A44: Dynamic Model: Estimating the Impact of GATT/WTO on PM10 by Regime Type
eststo clear
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto i.year i.region  if polity2 > 6 & polity2 != ., robust 
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != ., robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != ., robust

eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto i.year i.region  if polity2 <= 6 & polity2 !=., robust 
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != ., robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != ., robust

eststo: xi: xtreg log_PM10 L1log_PM10 L1wto i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust

eststo: xi: xtreg log_PM10 L1log_PM10 L1wto i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg log_PM10 L1log_PM10 L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
*esttab using "PM10Regime(Dynamic).tex", booktabs label title(Dynamic Model: Estimating the Impact of GATT/WTO on $PM10$ by Regime Type\label{PM10(Dynamic)Demo}) keep (L1log_PM10  L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1log_PM10 )  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A45: Difference in Difference: Estimating the Impact of GATT/WTO on Coal Energy Use by Regime Type
eststo clear
eststo: xi: xtreg electricity_coal_production gatt_wto i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != .,  fe robust

eststo: xi: xtreg electricity_coal_production gatt_wto i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_coal_production gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != .,  fe robust

eststo: xi: xtreg electricity_coal_production wto i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_coal_production wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_coal_production wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust

eststo: xi: xtreg electricity_coal_production wto i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_coal_production wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_coal_production wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
*esttab using "coalRegime(DID).tex", booktabs label title(Difference in Difference: Estimating the Impact of GATT/WTO on  Coal Energy Use by Regime Type \label{coal(DID)Demo}) keep(gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A46: Dynamic Model: Estimating the Impact of GATT/WTO on Coal Energy Use by Regime Type
eststo clear
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto i.year i.region  if polity2 > 6 & polity2 != ., robust 
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != ., robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != ., robust

eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto i.year i.region  if polity2 <= 6 & polity2 !=., robust 
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != ., robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != ., robust

eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust

eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1log_gdp_pc i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_coal_production L1electricity_coal_production L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
*esttab using "coalRegime(Dynamic).tex", booktabs label title(Dynamic Model: Estimating the Impact of GATT/WTO on Coal Energy Use by Regime Type\label{coal(Dynamic)Demo}) keep (L1electricity_coal_production  L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1electricity_coal_production )  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A47: Difference in Difference: Estimating the Impact of GATT/WTO on Renewable Energy Use by Regime Type
eststo clear
eststo: xi: xtreg electricity_production_renew gatt_wto i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != .,  fe robust

eststo: xi: xtreg electricity_production_renew gatt_wto i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != .,  fe robust
eststo: xi: xtreg electricity_production_renew gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != .,  fe robust

eststo: xi: xtreg electricity_production_renew wto i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_production_renew wto log_gdp_pc log_pop i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_production_renew wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 > 6 & polity2 != . & year>1994,  fe robust

eststo: xi: xtreg electricity_production_renew wto i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_production_renew wto log_gdp_pc log_pop i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
eststo: xi: xtreg electricity_production_renew wto log_gdp_pc log_pop life_expectancy urban_pop_per i.year if polity2 <= 6 & polity2 != . & year>1994,  fe robust
*esttab using "renewRegime(DID).tex", booktabs label title(Difference in Difference: Estimating the Impact of GATT/WTO on  Coal Energy Use by Regime Type \label{renew(DID)Demo}) keep(gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*"  ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps

*Table A48: Dynamic Model: Estimating the Impact of GATT/WTO on Renewable Energy Use by Regime Type
eststo clear
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto i.year i.region  if polity2 > 6 & polity2 != ., robust 
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != ., robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != ., robust

eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto i.year i.region  if polity2 <= 6 & polity2 !=., robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1log_gdp_pc L1log_pop i.year i.region if polity2 <= 6 & polity2 != ., robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1gatt_wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != ., robust

eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1log_gdp_pc L1log_pop i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 > 6 & polity2 != . & year>1994, robust

eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1log_gdp_pc L1log_pop  i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
eststo: xi: xtreg electricity_production_renew L1electricity_production_renew L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per i.year i.region if polity2 <= 6 & polity2 != . & year>1994, robust
*esttab using "renewRegime(Dynamic).tex", booktabs label title(Dynamic Model: Estimating the Impact of GATT/WTO on Coal Energy Use by Regime Type\label{renew(Dynamic)Demo}) keep (L1electricity_production_renew  L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per) order (L1gatt_wto L1wto L1log_gdp_pc L1log_pop L1life_expectancy L1urban_pop_per L1electricity_production_renew )  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Year FE=_Iyear_*" "Region FE=_Iregion_*" ) nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps


*Table A49: Probit: Estimating the Impact of GATT/WTO on FIT by Regime Type
eststo clear
eststo:  xtprobit FIT1 gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per time time2 time3 i.region_id if polity2 > 6 & polity2 != .
eststo:  xtprobit FIT1 gatt_wto log_gdp_pc log_pop urban_pop_per time time2 time3 i.region_id if polity2 <= 6 & polity2 != .
eststo:  xtprobit FIT1 wto log_gdp_pc log_pop life_expectancy urban_pop_per time time2 time3 i.region_id if polity2 > 6 & polity2 != .
eststo:  xtprobit FIT1 wto log_gdp_pc log_pop urban_pop_per time time2 time3 i.region_id if polity2 <= 6 & polity2 != .
*esttab using "FIT(probit)Regime.tex", booktabs label title(Probit: Estimating the Impact of GATT/WTO on FIT by Regime Type \label{FITDemo}) keep (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per)  mtitles ("" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Time Trend=time" )  nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps


*Table A50: Ordered Probit: Estimating the Impact of GATT/WTO on Emissions Standard by Regime Type
eststo clear
eststo: xi: xtoprobit auto_emissions_reg gatt_wto L1auto_emissions_reg time time2 time3 i.region if polity2 > 6 & polity2 != ., vce(robust)
eststo: xi: xtoprobit auto_emissions_reg gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if polity2 > 6 & polity2 != ., vce(robust)
eststo: xi: xtoprobit auto_emissions_reg gatt_wto L1auto_emissions_reg time time2 time3 i.region if polity2 <= 6 & polity2 != ., vce(robust)
eststo: xi: xtoprobit auto_emissions_reg gatt_wto log_gdp_pc log_pop life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if polity2 <= 6 & polity2 != ., vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto L1auto_emissions_reg time time2 time3 i.region if polity2 > 6 & polity2 != . & year>1994, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto log_gdp_pc log_pop life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if polity2 > 6 & polity2 != . & year>1994, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto L1auto_emissions_reg time time2 time3 i.region if polity2 <= 6 & polity2 != . & year>1994, vce(robust)
eststo: xi: xtoprobit auto_emissions_reg wto log_gdp_pc log_pop life_expectancy urban_pop_per L1auto_emissions_reg time time2 time3 i.region if polity2 <= 6 & polity2 != . & year>1994, vce(robust)
*esttab using "autoemissionRegime.tex", booktabs label title(RE Ordered Probit: Estimating the Impact of GATT/WTO on Auto Emissions Regulation by Regime Type \label{autoDemo}) keep(gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per L1auto_emissions_reg) order (gatt_wto wto log_gdp_pc log_pop life_expectancy urban_pop_per L1auto_emissions_reg)  mtitles ("" "" "" "" "" "" "" "" "" "" "" "")  stats(N_g N, labels(Countries Observations) fmt(0 0)) indicate("Time Trend=time" "Region FE=_Iregion_*") nodepvars se(3) b(3) replace star(* 0.10 ** 0.05 *** 0.01) compress nogaps





